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Abstract. We have performed shell-model calculations of the half-lives and neutron-branching probabilities 
of the r-process waiting point nuclei at the magic neutron number A'' = 82. These new calculations use 
a larger model space than previous shell model studies and an improved residual interaction which is 
adjusted to recent spectroscopic data around A = 130. Our shell-model results give a good account of 
all experimentally known half-lives and Q^-values for the N = 82 r-process waiting point nuclei. Our 
half-life predictions for the N = S2 nuclei with Z = 42-46 agree well with recent estimates based in the 
energy-density functional method. 
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1 Introduction 

The astrophysical r-process produces about half of the 
heavy elements in the Universe by a sequence of fast neutron- 
capture reactions interrupted by photodissociations and 
followed by (3 decays, running through extremely neutron- 
rich nuclei far off the valley of stability [HIS] . The /3 decays 
of the r-process waiting points, associated with nuclei with 
magic neutron numbers iV = 50, 82 and 126 play a cru- 
cial role for the r-process dynamics and elemental abun- 
dance distributions [3]. Despite their importance only a 
few half-lives of waiting points with magic neutron num- 
bers = 50 and 82 are known experimentally [4j[5j, while 
no experimental data exist yet for the N = 126 waiting 
points. This paper is concerned with the = 82 wait- 
ing points, and fortunately here the half-lives of ^^^In, 
^■^''Cd and ^^^Ag have been measured and serve as strin- 
gent constraints for models which have to be used to pre- 
dict the unknown half-lives of the other waiting points and 
r-process nuclei. 

Beta-decays are notoriously difficult to model as they 
are determined by the weak low-energy tails of the Gamow- 
Teller strength distribution, mediated by the operator <tt_ . 
There have been several previous estimates for the half- 
lives of the = 82 waiting points based on the Quasi- 
particle Random Phase Approximation on top of semi- 
empirical global models [6l[7] , the energy-density functional 
(DF3) method [S] or the Hartree-Fock-Bogohubov (HFB) 
method [Qj. A comparison to the data show that the pre- 
dicted half-lives are in the right order of magnitude, but 
are often somewhat too long. This might imply that the 
models underestimate the correlations among nucleons whichimproved study. 



pull down the Gamow- Teller (GT) strength to low exci- 
tation energies. It is well known that the interacting shell 
model is the method of choice to describe the Gamow- 
Teller distribution in nuclei [IO j ll H [T2] and in fact, the best 
agreement with the measured half-lives has been achieved 
within a shell model approach |13j . 



This shell model calculation as well as the other the- 
oretical approaches have been challenged recently by the 
experimental determination of the excitation energy of the 
first 1+ state in ^■^"in, which carries most of the GT_ 
strength of the ^^°Cd decay [13]. Here, gamma rays ob- 
served in the beta decay of ^^°Cd [14] yield an excita- 
tion energy of = 2.16 MeV, in contrast to the shell 
model predictions which place this state at noticeably 
lower energy, = 1.4 MeV [M] and 1.5 MeV [13:. As 
the halflife has a strong energy dependence which approx- 
imately scales like {Qi3 — Ex)^, where Qp = Mi — Mf is the 
difference of the masses of the parent and daughter nuclei, 
respectively, the misplacement of the 1"*" excitation energy 
has been fortuitiously cancelled in the shell model calcu- 
lation 113] by the use of too small a Q/j value which was 
taken from the Duflo-Zuker mass model yielding 7.56 MeV 
[15], while the experimental Qp value is 8.34 MeV. Thus 
the shell model calculation of [13] despite its successful 
desription of the halflives of the N = 82 waiting point 
nuclei, has to be improved to account also for the recent 
experimental structure data concerning the respective de- 
cays. It is the aim of this manuscript to present such an 
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We have performed large-scale shell model calculations us- 
ing the code ANTOINE [16 . We have improved the pre- 
vious shell model study of the half-lives of the iV = 82 
r-process nuclei |13| in several ways. At first, we used a 
larger model space which includes the 037/2, 1^3/2,5/2: 2si/2, 
orbitals outside the = 40 core for neutrons, thus as- 
suming a closed A^ = 82 shell configuration in the parent 
nucleus. For protons our model space was spanned by the 
059/2,7/2, l'^3/2. 5/2: 2si/2 orbitals. Thus our model space 
avoids spurious center-of-mass excitations by omitting the 
hii/2 orbit for protons and the 59/2 orbit for neutrons. 
Secondly, we have performed the calculations of the par- 
ent ground states and the GT strength distributions in 
the daughter nucleus for the A^ = 82 parents with charge 
numbers Z = 43 — 49 including all possible correlations 
within the defined model space, a clear improvement with 
respect to the calculations in ref. [13]. But most impor- 
tantly we have modified the residual interaction adopted 
in our studies to reproduce relevant experimental nuclear 
structure information. These modifications were guided 
by the observation that for the nuclei of interest here, 
the halflives are dominated by Gamow- Teller transitions 
to states at low excitation energy in the daughter nuclei 
and that these transitions are mainly determined by a sin- 
gle transition matrix element in which a 37/2 neutron is 
changed into a gg/2 proton. This implies that the transi- 
tion matrix elements are quite insensitive to the relative 
position of the 37/2 neutron orbit but its position deter- 
mines the excitation energy of the daughter states and 
the (5/3 value for the evaluation of the half-lives. Start- 
ing from the interaction used in the previous shell-model 
calculations [13j we have implemented several monopole 
modifications aiming to reproduce the known Qp values 
of i3iln [17] and "°Cd [M], the position of the /in/2 neu- 
tron orbit [17] and the experimental excitation energy of 
the first 1+ state in ^"^"in (for details see ref. [H]). This 
has resulted in two effective interactions: one that allows 
proton excitations from the pi/2 orbit (that means uses 
a ^^Sr core) and another one where these excitations are 
supressed (uses a ^°Zr core). Fig. [1] compares our calcu- 
lated low-energy spectrum of ^^°In with the data and find 
good agreement for the excitation energies of the first 3"'' 
state. Importantly our improved shell model calculations 
also reproduces the unexpectedly high excitation energy 
of the first 1+ state, which is of key importance for the cal- 
culation of the ^^°Cd halflife. However, the low energy 5"'' 
state and the states with possible assigments 0^ and 1~ 
are missed by the calculation that uses a ^"Zr core. These 
are reproduced by the interactions that includes the pi ji 
orbital which is energetically relatively close to the 179/2 
orbital, but carries the opposite parity. Once, the pi/2 or- 
bital is included in the model space we are forced to trun- 
cate the number of protons excited across the 39/2 shell 
gap. The calculations shown at the right of figure [1] where 
performed allowing for 6 protons to be excited from the 
Pi/2 and (79/2 orbitals. Shell model calculations performed 
within the same model space and using the same inter- 
action also reproduce the recently measured spectrum of 
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Fig. 1. Comparison of low-energy shell model spectrum for 
^^°In with the data. The shell model spectrum of the left is 
calculated for the model space assuming a ^°Zr core, the one 
on the right for the model space with a **Sr core. Details of 
the calculations are given in the text. 



^^°Cd [IHl. It is particularly noteworthy that this includes 
the excitation energy of the first excited 2+ state which 
is calculated at = 1.325 MeV in close agreement with 
the experimental value of 1.346 MeV, hence not confirm- 
ing the tentative suggestion that this state would reside at 
Ex — 0.957 MeV which had been interpreted as an onset 
of shell quenching already in ^^°Cd [50]. 

As stated above, our calculation of the halflives are 
based on the valence space on top of the ^°Zr as this al- 
lows for untruncated calculations and the changes in the 
computed values of the half-lives are negligible if we adopt 
the model space with the '^^Sr core (and the appropriate 
interaction). 



3 Results 

As discussed above, /3 half-lives depend sensitively on the 
value. Unfortunately this important quantity is ex- 
perimentally not known for most of the relevant N = 82 
r-process nuclei and hence has to be estimated based on 
theoretical models or by systematic extrapolations from 
data for neighboring nuclei. Our Qp values were calcu- 
lated from the energies of the isobaric analog states with 
a systematic correction for the Coulomb displacement en- 
ergies [23]. In Fig. [2] we compare the shell model (5/3-values 
with experiment, the Audi-Wapstra systematics [21j and 
other theoretical calculations. As stated above, our shell 
model calculation is tuned by proper monopole modifica- 
tions to reproduce the experimental Q/s value for ^■^"Cd. 
However, we find also a good agreement to the other Qp 
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Fig. 2. Comparison of shell model Qp values for the N — S2 
isotones to the Audi-Wapstra systematics [21] and to predic- 
tions of other models: FRDM [B], ETFSI-Q [22], DF3 [8] and 
Duflo-Zuker |15|. 



values as given in the compilation of Audi-Wapstra [IJ. 
Thus as the shell model Qp values are close to the exper- 
imental values and agree with those of most other models 
within the theoretical uncertainties, we will in the follow- 
ing adopt the shell model Q^g values in our calculation of 
the half-lives. However, the exception are the predictions 
of the quenched Extended Thomas-Fermi with Strutinski 
Integral (ETFSI-Q) model [22] which predicts Q/j-values 
which are noticeably larger than those of the other mod- 
els, most noticeably for the most neutron deficient = 82 
isotones. 

The final ingredients of our half- life calculations are the 
GT_ strength functions. We have calculated them within 
our shell model approach using the Lanczos method with 
60 iterations for each possible final J-value. The calcu- 
lated GT_ strength functions for the A^ = 82 nuclei from 
124Mo {Z = 42) to i3iln {Z = 49) are shown in Figs.[5][nj 
For even-even nuclei the GT transitions lead to JJ = 1+ 
states in the daughters, while for the odd-A nuclei the fi- 
nal states can have Jf = Ji — 1, Ji, Ji + 1, where Ji is the 
angular momentum of the parent ground state. For the 
calculation of the half-lives, the choice of the appropriate 
Gamow- Teller quenching factor is an important issue. It 
is wellknown that shell model calculations reproduce the 
GT strength distributions (total strength and fragmen- 
tation) very well within complete OHlo calculations (i.e. 
model spaces which include a complete major oscillator 
shell) , if the GT operator is quenched by a constant fac- 
tor. This factor has been determined for sd shell ;24;„.10j, 
where it is 0.77, and the pf shell, where it is 0.74 pS] . 
Unfortunately the appropriate constant for nuclei in the 
mass A — 130 region has not yet been determined. Thus 
we adopt a quenching factor of 0.71 adjusted to reproduce 
the experimental half-life of ^^°Cd {ti^2 = 162 ± 7 ms). 

Here one word about first-forbidden transitions is in 
order. While a study by MoUer and collaborators, based 
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Fig. 3. Comparison of half-lives of the N — 82 isotones as 
calculated in the FRDM ^6^, HFB [5], and the present shell 
model approaches with data [4lll4j. 



on the FRDM/QRPA for GT transitions and the sta- 
tistical gross theory for forbidden theories, indicates siz- 
able contributions of first-forbidden transitions to the half- 
lives, the only consistent microscopic treatment of GT and 
first-forbidden transitions for the A^ = 82 half-lives by 
Borzov, based on the energy density-functional method, 
implies that forbidden transitions accelerate the halflives 
only slightly by about 10% or less. Based on the later re- 
sult our halflives calculated purely from GT transitions are 
meaningful, as a small, but roughly constant contribution 
of first-forbidden transition can be absorbed into a modi- 
fication of the quenching factor which would be 0.73 (very 
close to the standard quenching factor), if first-forbidden 
transitions yield a 10% contribution. 

Our shell model half-life for ^^^In (260 ms) agrees also 
with the measured values (280 ±30 s). However, we over- 
estimate the one for ^^^Ag slightly (70 ms to be compared 
with 461^9 ^) [26,4J. The present shell- model half-lives are 
longer for nuclei with Z < 47 than the ones computed 
in reference |13j (see table [T]) . There are two reasons for 
the change in half-lives. First, for nuclei with Z = 42-44 
the Q^-values obtained in the present shell- model calcula- 
tions are around 1 MeV smaller than the Duflo-Zuker val- 
ues used in ref. [T^ (see fig. dj). Second, for all the nuclei 
the low lying Gamow- Teller strength is shifted upwards 
around 0.5 MeV in excitation energy. This shift results in 
longer half-lives even if the Q/j-value were unmodified as 
it is the case in ^^^Ag. This shift is due to the monopole 
modifications introduced in the interaction used for the 
present calculations to increase the excitation energy of 
the 1"*" state in ^■^''in from 1.6 MeV [T3| to the experimen- 
tal value of 2.12 MeV [H]. 

In Fig. [3] and Table [T] our results are compared with 
those of other theoretical models and the data. The cur- 
rent half-lives are quite similar to those obtained within 
the HFB model of |£j and are also in good agreement 
with the DF3-QRPA approach [57]. The latter is reassur- 
ing as the DF3-QRPA model also reproduces the half-lives 
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Table 1. Comparison of the present shell model half- lives and 
the ones of reference [13) with experiment [26ll4l[T4]. All half- 
lives are in ms. 



Nucleus 




Half-Life (ms) 


Expt. 


Theor. 


Theor. (ref [13]) 


i^iln 


280 ± 30 


260 


177 


"OCd 


162 ±7 


162 


146 


129 Ag 




70 


35.1 


128pd 




46 


27.3 


i^^Rh 




27.65 


11.8 


i^f^Ru 
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9.6 
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Fig. 12. Comparison of the shell model neutron separation 
energies for the = 81 isotones to the data [21] and predic- 
tions of other models: FRDM [6], Duflo-Zuker [15] and ETFSI- 
Q i22,. 



of other spherical nuclei in the vicinity of iV = 82 quite 
well [28 . Our half-lives show a mild odd-even effect, with 
the half-lives of the even-even waiting point nuclei slighly 
enlarged with respect to the neighboring odd-A nuclei. 
This is caused by a partial cancellation of the odd-even 
staggering in the Q/j values (Fig. [J) by the larger excita- 
tion energies of the lowest GT states in the odd-^ daugh- 
ter nuclei (Figs. IMTTjl. 

As shown in Fig. [2] the various models predict Qjs val- 
ues within a variation of about 1 MeV, which, due to the 
strong energy dependence of the phase space, translates 
into the largest uncertainties of the (3 half-lives. To esti- 
mate this uncertainty we have recalculated the half-lives 
by replacing the shell model values by those of the 
different models and find half-lives which agree with the 
present ones within a factor of two or better, except for 
the larger ETFSI-Q Qp values which result in significantly 
smaller half-lives for the most proton-deficient nuclei. It is 
also worth noting that by replacing our shell model Qp 
value for ^^^Ag by the systematic Audi-Wapstra value, we 
find agreement with the experimental half-life within the 
uncertainties of the systematic Qfj value. 

As the neutron separation energies in the daughter nu- 
clei are quite small, some of the GT_ strengths resides 
actually at energies above the neutron emission thresh- 
old and (3 decays to these states are followed by neutron 
emission. The probability for /3-delayed neutron emission 
depends sensitively on the neutron separation energies Sn 
which are not known experimentally for most of the nuclei 
of interest here and have to be estimated theoretically. In 
Fig. [12] we compare our shell model neutron separation en- 
ergies for the TV = 81 daughter nuclei with those of various 
models. Importantly for ^^^Sn and ^■^"in the neutron sepa- 
ration energies are known experimentally and our present 
results agree quite nicely with the data. The neutron sepa- 
ration energies predicted in the FRDM model [6] are larger 
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than all other model predictions and also exceed the ex- 
perimental values for ^^^Sn and ^^°In. Our shell model re- 
sults predict a quite similar slope of the Sn values as found 
in the Duflo-Zuker mass model [TS] yielding quite sizable 
neutron separation energies even in the proton-deficient 
nuclei ^^^Ru and ^^'^Tc. This is in difference to the ETFSl- 
Q values [22] which predict a pronounced weakening of the 
Sn values with increasing neutron excess. A quite similar 
behavior is found if one compares the 2-neutron separation 
energies of the = 82 r-process nuclei (Fig. [13]). Again, 
our shell model results agree with the available data (^'^^In 
and "'^^°Cd) and show a significantly slower decrease of the 
S2n values than predicted by the ETFSl model, while the 
FRDM and Duflo-Zuker models yield 2-neutron separa- 
tion energies which agree reasonably well with the shell 
model ones. 



Using our shell model GT_ strength functions and Sn 
values we have calculated the probability P„ that the /? 
decay is accompanied by the emission of (at least) one neu- 
tron, defined as the relative probability of the /3-decay rate 
above the neutron emission threshold Sn- The results are 
shown in Fig. [TJ] indicating that within our shell model 
study most of the (3 decays go to states below the neutron 
threshold for the even-even parent nuclei, yielding proba- 
bilities for /3-delayed neutron emission of 40% or less. The 
shell model predicts quite a strong odd-even staggering in 
the Pn values, indicating P„ values of 80% or larger for 
i27Rh and i25tc. The FRDM model |6] also predicts an 
odd-even dependence in the neutron emission probabili- 
ties, however, this effect is somewhat smaller than for the 
shell model values. For ^'^°Cd, the P„ value is known ex- 
perimentally [14j and it agrees with the FRDM and shell 
model predictions. 



4 Conclusion 

We have recalculated shell model half-lives and neutron 
emission probabilities for the A^ = 82 waiting point nuclei 
in the r-process, improving a previous shell model study by 
enlargement of the model space and by modification of the 
residual interaction which reproduces recent spectroscopic 
findings for nuclei in this regime of the nuclear chart. In 
particular our modified calculation reproduces the unex- 
pectedly high excitation energy of the first 1+ state in 
^^°In |14j. This, as a good description of the Qp value, is 
crucial as the GT transition to this low-lying state domi- 
nates the ^'^"Cd half-life. We find good agreement with the 
experimentally known half-lives for the Z = 48,49 nuclei 
and overestimate the one of ^^^Ag slightly. 

As in most other models [5J[ni[7] the shell model half- 
lives have been computed solely on the assumption of 
an allowed Gamow- Teller transition. The importance of 
forbidden transitions for the half-lives of the A^ = 82 
waiting point nuclei is somewhat controversial. Using the 
gross theory for the forbidden transitions, MoUer [2^ finds 
non-negligible contributions, while Borzov, adopting the 
QRPA approach within the Fermi-liquid theory, concludes 
that forbidden transitions play only a minor role [27]. In 
view of these differences, shell model calculations of the 
forbidden /3 decays of the A^ = 82 r-process waiting point 
nuclei are desirable. 

I. N. Borzov is supported by a grant from the German DFG 
under contract number 436 RUS 113/907/0-1. 
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